Spectral Analysis of Virus Spreading in Random Geometric Networks 



Victor M. Preciado and Ali Jadbabaie 



Abstract — In this paper, we study the dynamics of a viral 
spreading process in random geometric graphs (RGG). The 
spreading of the viral process we consider in this paper is 
closely related with the eigenvalues of the adjacency matrix 
of the graph. We deduce new explicit expressions for all the 
moments of the eigenvalue distribution of the adjacency matrix 
as a function of the spatial density of nodes and the radius of 
connection. We apply these expressions to study the behavior of 
the viral infection in an RGG. Based on our results, we deduce 
an analytical condition that can be used to design RGG's in 
order to tame an initial viral infection. Numerical simulations 
are in accordance with our analytical predictions. 

I. Introduction 

The analysis of spreading processes in large-scale complex 
networks is a fundamental dynamical problem in network 
science. The relationship between the dynamics of epi- 
demic/information spreading and the structure of the under- 
lying network is crucial in many practical cases, such as 
the spreading of worms in a computer network, viruses in a 
human population, or rumors in a social network. Several 
papers approached different facets of the virus spreading 
problem. A rigorous analysis of epidemic spreading in a 
finite one-dimensional linear network was developed by 
Durrett and Liu in [3]. In [10], Wang et al. derived a sufficient 
condition to tame an epidemic outbreak in terms of the 
spectral radius of the adjacency matrix of the underlying 
graph. Similar results were derived by Ganesh et al. in [4], 
establishing a connection between the behavior of a viral 
infection and the eigenvalues of the adjacency matrix of the 
network. 

In this paper, we study the dynamics of a viral spreading 
in an important type of proximity networks called Random 
Geometric Graphs (RGG). RGG's consist of a set of vertices 
randomly distributed in a given spatial region with edges 
connecting pairs of nodes that are within a given distance 
r from each other (also called connectivity radius). In this 
paper, we derive new explicit expressions for the expected 
spectral moments of the random adjacency matrix associated 
to an RGG. Our results allow us to derive analytical condi- 
tions under which an RGG is well-suited to tame an infection 
in the network. 

The paper is structured as follows. In Section II, we 
describe random geometric graphs and introduce several 
useful results concerning their structural properties. We also 
present the spreading model in [10] and review an important 
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result that relates the behavior of an initial infection with the 
spectral radius of the adjacency matrix. In Section III, we 
study the eigenvalue spectrum of random geometric graphs. 
We derive explicit expressions for the expected spectral 
moments in the case of one- and two-dimensional RGG's. 
In Section IV, we use these expressions to study the spectral 
radius of RGG's. Our results allow us to design RGG's 
with the objective of taming epidemic outbreaks. Numerical 
simulations in Section IV validate our results. 

II. Virus Spreading in Random Geometric Graphs 

In this section, we briefly describe random geometric 
graphs and introduce several useful results concerning their 
structural properties (see [7] for a thorough treatment). We 
then describe the spreading model introduced in [10] and 
show how to study the behavior of an infection in the network 
from the point of view of the adjacency eigenvalues. 

A. Random Geometric Graphs 

Consider a set of « nodes, V„ — {vi,...,v n }, respectively 
located at random positions, % n = {xi,...,x„}, where X; 
are i.i.d. random vectors uniformly distributed on the d- 
dimensional unit torus, T d . We use the torus for convenience, 
to avoid boundary effects. We then connect two nodes v,, v; <G 
V„ if and only if x,-— xy|| < r, where r is the so-called 
connectivity radius. In other words, a link exists between 
v, and Vj if and only if Vy lies inside the sphere of radius 
r(n) centered at v,-. We denote this spherical region by S;(r), 
and the resulting random geometric graph by G(% n ;r). We 
define a walk of length k from vo to as an ordered set of 
(possibly repeated) vertices (vo, Vi , v#) such that V; ~ v, + i , 
for i = 0, 1, ...,k— 1; if = Vo the walk is said to be closed. 

The degree d\ of a node V; is the number of edges 
connected to it. In our case, the degrees are identical random 
variables with expectation [7]: 

E[di] = «V (c V, (1) 

where VW is the volume of a <f -dimensional unit sphere, 
v (d) = n d/^ /Y(d/2 + \), and T(-) is the Gamma function. 
The clustering coefficient is a measure of the number of 
triangles in a given graph, where a triangle is defined by 
the set of edges {(i,j) , (j,k) , (k,i)} such that i ~ j ~ k^ i. 
For one- and two-dimensional RGG's we can derive an 
explicit expression for the expected number of triangles, 
E[f,-], touching a particular node v, (details are provided in 
Section III). 

The adjacency matrix of an undirected graph G, denoted 
by A{G) = [ay], is defined entry-wise by ay = 1 if nodes 
i and j are connected, and ay = otherwise. (Note that 



a,, = for simple graphs.) Denote the eigenvalues of a n x n 
symmetric adjacency matrix A(G) by X\ < ... < X n . The k-th 
order moment of the eigenvalue spectrum of A(G) is defined 
as: 

JV 1=1 

(which is also called the k-t\\ order spectral moment). 

We are interested in studying asymptotic properties of the 
sequence G(% n \r («)) for some sequence {r(n) : n £ N}. In 
[7], two particularly interesting regimes are introduced: the 
thermodynamic limit with nr(nf — > a E (0,°°), so that the 
expected degree of a vertex tends to a constant, and the 

connectivity regime with r (n) — > yl — |^ ) with a constant 
y, so that the expected degree of the nodes grows as clogn. 
In this paper, we focus on studying the spectral moments 
in the connectivity regime. In Section III, we derive explicit 
expressions for the expected spectral moments of G(Xn',r„) 
for any network size «. We then use this information to bound 
the spectral radius of the adjacency matrix of G{% n ;r{n)). 

B. Spectral Analysis of Virus Spreading 

In this section, we briefly review an automaton model 
that describes the dynamics of a viral infection in a spe- 
cific network of interactions. This model was proposed and 
analyzed in [10], where a connection between the growth of 
an initial infection in the network and the spectral radius of 
the adjacency matrix was established. This model involves 
several parameters. First, the infection rate j3 represents the 
probability of a virus at an infected node i spreading to 
another neighboring node j during a time step. Also, we 
denote by 8 the probability of recovery of any infected node 
at each time step. For simplicity, we consider j3 and 5 to be 
constants for all the nodes in G. We also denote by pi [k] the 
probability that node ; is infected at time k. The evolution 
of the probability of infection is modeled by means of the 
following system of non-linear difference equation: 

Pi [k+ 1] = [1 - n (1 -p Pj [k])} + (1 - 8) p t [k] , (2) 

for i= 1, ...,«, where jYi denotes the set of nodes connected 
to node i. We are interested in studying the dynamics of the 
system for a low-density level of infection, i.e., j3 pj [k] <IC 
1. In this regime, a sufficient condition for a small initial 
infection to die out is [10]: 

A max (A(G))<-. (3) 

One can prove that |3) is a sufficient condition for local 
stability around the disease-free state. Thus, we can use 
condition (O to design networks with the objective of taming 
initial low-density infections. 

III. Spectral Analysis of Random Geometric 
Graphs 

In this paper, we study the eigenvalue distribution of the 
random adjacency matrix associated to G{% n ;r (n)) for n — > 
oo. In this section, we characterize eigenvalue distribution 



using its sequence of spectral moments. In our derivations, 
we use an interesting graph-theoretical interpretation of the 
spectral moments [1]: the k-th spectral moment of G is 
proportional to the number of closed walks of length k in G. 
This result allows us to transform the algebraic problem of 
computing spectral moments of the adjacency matrix into the 
combinatorial problem of counting closed walks in the graph. 
In the following subsection, we compute the expected value 
of the number of closed walks of length k in G(Xn\ r {n)). 

A. Spectral Moments of One-Dimensional RGG 's 

As we mentioned above, we can compute the £-th spectral 
moment of a graph by counting the number of closed walks 
of length k. In the case of an RGG G{% n ;r{ri)), this number 
is a random variable. In this subsection, we introduce a novel 
technique to compute the expected number of closed walks 
of length k. For clarity, we introduce our technique for the 
first three expected spectral moments k= 1,2,3. We then use 
these results to induce a general expression for higher-order 
moments in one-dimensional RGG's. 

The first-order spectral moment is equal to the number 
of closed walks of length k = 1. Since G(^„;r) is a simple 
graphs with no self-loops, we have that m\{G(% n \r)) is a 
deterministic quantity equal to 0. 

We now study the expected second moment, 
E [ni2 {G{x„\r))\, by counting the number of closed 
walks of length two. In simple graphs, the only possible 
closed walks of length two are those that start at a given 
node v,, visit a neighboring node vj E jfi, and return back 
to V/. Hence, the number of closed walks of length two 
starting at v,- is equal to d,. Thus, from ((TJ, we have 

E\m 2 } = -YMdi\=nV [d) r d , 

where this result is valid for any dimension d > 1 . 

The third spectral moment is proportional to the number 
of closed walks of length three in the graph. We now derive 
an expression for the expected number of triangular walks 
starting at a given node v,- in a one-dimensional RGG. Since 
all nodes are statistically equivalent, our result is valid for 
any other starting node. For simplicity in our calculations, 
we consider that v, is located at the origin. A triangular walk 
starting at node v; exists if and only if there exist two nodes 
vy anv Vh such that |xy| < r, |jc>fc| < r, and pcjt — Xj\ < r. Also, 
since the random distribution of vertices on T 1 is uniform 
(with density «), the probability of nodes Vj and v* being 
respectively located in the differential lengths [xj + dxj) and 
[x/c + dxk) is equal to n 2 dxjdx^. Hence, one can compute the 
expected number of triangular walks starting at node v, as 

/ / n 2 dxjdxk, 

J J(x Jr x k )eH 2 (r( n )) 

where 

H 2 (r) = {(xj,x k )eT 2 s.t. \xj\<r, (4) 
I**-*/! < r,\x k \ < r} . 



Thus, E[f,-] can be computed as n 2 Vol [H2 (r(n))] (where 
Vol(H) denotes the volume contained by the polyhedron H.) 
Notice that H2 (r) can be defined by a set of linear inequalies; 
hence, H2 (r) is a convex polyhedron that depends on r. 
Furthermore, the set of linear inequalities in dU presents a 
homogeneous dependency with respect to the parameter r. 
Therefore, we can write Vol (H2 (r)) as r 2 Vol (H2 (1)). Finally, 
one can easily compute the volume of H2 (1) to be equal 
to 3. Thus, the expected third spectral moment of a one- 
dimensional RGG is given by 



E [m 3 ] 



11 «™ 



3n 2 r 2 . 



i=i 



In the following, we extend the above technique to com- 
pute higher-order expected spectral moments. Denote by 

(k) 

W t the number of closed walks of length k starting at node 

we derive the following 



(k) 



vi in G{% n ;r{n)). Regarding W t 
result. 

Theorem 1: The expected number of closed walks of 
length k, W- , in a random geometric graph, G(% n \r), on 
T 1 is given by 

" /(*) 



E 



w.y 



(nr) 



2(k- 



k-2 



E k - 



where E^- 1 ; are the Eulerian numbers Q. 

Proof: Consider a particular closed walk, w* = 
(vi,V2,V3,...,v^,vi), of length k starting and ending at node 
vi (which we locate at zero for computational convenience). 
A walk Wjt exists if and only if there exists a set of k— 1 
nodes, {v2, V3, vj} , such that |jci| < r, \xj+\—Xj\ < r for 
j = 2, ...,k— 1, and \x^\ < r. Since the distribution of vertices 
on T 1 is uniform (with density n) one can compute the 
expectation of Wj as 

" ,(*) 



w 



(x 2 ,...,x k )eH k _i(r{n)) 



n k 1 dx2---dxk, 



where 



Hk-i(r) = {( 



V2,V3,- 

| x ;'+i ~ x j 
\x k \<r}. 



) e T*- 



s.t. 



< r for j = 2, 



|v 2 |<r, 

;k-l, 



(5) 



Thus, E[W^] can be computed as n k ~ l No\[Hi i _\ (r)], where 
Hk-i ( r ) is a convex polyhedron defined by a set of linear 
inequalities. Finally, note that the homogeneous structure of 
the system of linear inequalities defining H^-i (r) allows us 
to write Vol(/4_! (r)) = ^ Vol(H t _i (1)). Therefore, 



E 



(nr)*" 1 Vol (flin(l)). 



(6) 



The volume of (1) is a particular number, independent 
of the RGG parameters, i.e., n and r. Furthermore, we have 
found an explicit analytical expression for the volume of 
Hk (1) for any k > 1. Although we do not provide details of 

'The Eulerian number E(n,k) gives the number of permutations of 
{1,2, ...,n} having k permutation ascents [5]. 



our derivation, due to space limitations, an explicit expres- 
sion for the volume of Hk{\) is given by [9]: 

k- 

1 Ekj, (7) 



Voi(^(i)) = ^£ 



7=1 



k 

7-1 



where E^i denotes the Eulerian numbers. Substituting (0 in 
(O we obtain the statement of our lemma. ■ 
In [6], Lasserre proposed an algorithm to compute the 
volume of a polyhedron defined by a set of linear inequali- 
ties. We can use this algorithm to verify the validity of (O. 
Applying this algorithm to the set of inequalities in (0, we 
compute the following volumes for k = 1, 10: 

Hi =2, H 2 = 3, H 3 = 5.333..., H 4 = 9.58333..., 
H 5 = 17.6000..., H 6 = 32.70555..., H 7 = 61.3587..., 
H s = 115.947..., H 9 = 220.3238..., ffio = 420.825... 

These numerical values match perfectly with our analytical 
expression in Theorem Q] 

If nr(n) = i2(logn) (i.e., the average degree grows 
as log 71, or faster), one can prove that E^] = 
(1 +0(log- 1 «)) E[wj- k) ]. Hence, from © and Q, we 
have the following closed-form expression for the asymptotic 
expected spectral moments: 



E[m,]x(«r) 



1 



2(k- 



k-2 



k-i 
7-1 



E k - 



(8) 



In the following table we compare the analytical result 
in (HI) with numerical realizations of the empirical spectral 
moments. In our simulations, we distribute « = 1000 nodes 
uniformly in T 1 and choose a connectivity radius r = 0.01 
(which results in an average degree E[c/,] = 20). The second, 
third, and forth column in the following table represent 
the analytical expectations of the spectral moments, the 
empirical average of the spectral moments from 10 random 
realizations of the RGG, and the corresponding empirical 
typical deviation, respectively. 



k 


EH 


Empirical Average 


Typical Deviation 


1 





1.38e-16 


1.3e-15 


2 


20 


19.9326 


0.0976 


3 


300 


297.284 


4.3598 


4 


5,733 


5,956.30 


196.94 



Our numerical results present an excellent match with our 
analytical predictions. 

B. Spectral Moments of Two -Dimensional RGG's 

In this subsection, we derive expressions for the first three 
expected spectral moments of G(x„;r(«)) when the nodes 
are uniformly distributed in T 2 . The expressions for the 
first and second expected spectral moments are mi =0 and 
E['«2] = Knr 2 . The third spectral moment is proportional to 
the number of closed walks of length three in the graph. In 
the two-dimensional case, we count the number of triangular 
walks using a technique that we illustrate in Fig. 1. In this 
figure, we plot two nodes v, and Vj. The parameters p and 
§ in Fig. 1 denote the distance and angle between these 



Fig. 1. This figure illustrates the technique proposed in Section III.B to 
count the number of triangular walks in a two-dimensional RGG. 



two nodes, i.e., p = \\Xj— X;|| and = /(x 7 - — x,-). An edge 
between v, and vj exists if an only if Vj is located inside the 
circle 5,-(r). In this setting, the probability of existence of a 
triangle touching both v, and Vj is equal to the probability of 
a third node being in the shaded area A; (see Fig. 1). This 
area is the result of intersecting the circles 5, (r) and Sj(r), 
and the resulting probability is equal to n A/. The intersecting 
region A/ is a symmetric lens which area can be computed 
as a function of p and r as follows: 

„t \ I 2r 2 cos~ 



o, 



for p < r, 
for p > r. 



(9) 

Therefore, we can compute the expected number of triangles 
by integrating over the set of all possible positions of v/, i.e., 
77 e [0,r] and <j) g [0,2^), as follows 

E[n]= [' n 2 A,(p;r) p dp d<j). (10) 

Jp=0 J0=O 

After substituting (0 in ( flOt . we can explicitly solve the 
resulting integral to be 

3V3 S 



E u 



Ti- 



ll (nr 



:5.78(«r 2 ) 2 . 



(11) 



Consequently, we have the following expression for the third 
expected spectral moment E [/M3] = j Yli=\ E [f;] = E [?,-]. 
In the following, we extend the technique introduced above 

(k) 

to compute closed walks of arbitrary length. Denote by W/ 
the number of closed walks of length k starting at node v\ 
in G(% n ;r(n)). The idea behind our technique is illustrated 
in Fig. 2, where we represent a particular closed walk of 
length 6. We denote this walk by Wj = (vi,V2, vj_i , v^, vi). 
We define the following set of relative distances and angles 
between every pair of connected vertices: p, = ||x,-+i — x,-|| 
and (j>j = Z(x, + i — x,) for i= l,...,k — 2. We also define the 
following parameter 



k-2 
E Pi 6 



(12) 




Fig. 2. This figure illustrates the technique proposed in Section III.B to 
count the number of closed walks of length k in a two-dimensional RGG. 



(i = \J — 1) which is the resulting distance between nodes 
Vfc_i and vi given a particular set of distances and angles 
{('7j0i)};=l k-2 ( see Pig' ^)- In this setting, the conditional 
probability of existence of a walk = (vi , v%, Vfc_ 1 , v^, Vi) 
given the set of relative positions, {(n, 0i)}/=i i-2> I s e 1 ua l 
to the probability of being in the shaded area A/ in Fig. 
2. We have an expression for this area in (0, where p is 
defined in (fl2l) . Finally, we can compute the expectation of 

(k) 

W t by performing an integration over the set of all possible 
positions (i.e., pj € [0,r] and E [0,2n) for 7 =2,..., k— 1), 
as follows 



E 



W; 



(k) 



■/ 1 



k-l 

M{P\r) ]~[t7; dri d<p, 

J=2 



where J] = (p 2 , ...,Pk-i), 9 = (02,-,$fc-i), and C k _ 2 = 
{(rj,(p) : 77G[0,r] A '~ 2 and <p G [O^)*" 2 }. Although a 
closed-form for the above expression can only be computed 
for k < 3, we can always find a good approximation via 
numerical integration. For example, the integration for k = 4 
gives us E[W, (4) ] w 14.2511 (nr 2 ) 3 . 

In the following table, we compare our analytical results 
with numerical realizations of the empirical spectral mo- 
ments of a two-dimensional RGG. In our simulations, we 
distribute n = 1000 nodes uniformly on T 2 and choose a 
connectivity radius r = y^SO/nn rj 0.1784 (which results in 
an average degree E [*/,-] = 50). The second, third, and forth 
columns in the following table represent the analytical expec- 
tation of the spectral moments, the empirical average from 10 
random realizations, and the corresponding empirical typical 



deviation, respectively. 



k 


EH 


Empirical Average 


Typical Deviation 


1 





-9.2e-16 


l.le-15 


2 


50 


50.0820 


0.3908 


3 


1,464.1 


1,475.8 


37.3777 


4 


59,452 


60,127 


2,955.3 



Our numerical results present an excellent match with our 
analytical predictions. 

In the following section, we use the results introduced in 
this section to study the spreading of an infection in a random 
geometric network. 

IV. Spectral Analysis of Virus Spreading 

In this section, we use the expressions for the expected 
spectral moments to design random geometric networks to 
tame an initial viral infection in the network. In our design 
problem, we consider that the size of the network n and 
the parameters in (0, i.e., /3 and 5, are given. Hence, our 
design problem is reduced to studying the range of values of 
r for which the RGG is well-suited to tame an initial viral 
infection. 

A sufficient condition for local stability around the 
disease-free state was given in ([3]). Thus, we have to find 
the range of values of r for which the associated spectral 
radius Am ax is smaller than the ratio 5/j3. In the following 
subsection, we show how to derive an analytical upper 
bound for the spectral radius based on the expected spectral 
moments. 

A. Analytical Upper Bound for the Spectral Radius 

In order to upper-bound the spectral radius, we use 
Wigner's high-order moment method [11]. This method 
provides a probabilistic upper bound based on the asymptotic 
behavior of the A:-th expected spectral moments for large k. 
We present the details for a one-dimensional RGG, although 
the same technique can be applied to RGG's in higher 
dimensions. For a one-dimensional RGG in the connectivity 
regime, we derived an explicit expression for the expected 
spectral moments in ([8]). A logarithmic plot of Vol(/4(l)) 
for k = 1,2,..., 9 unveils that Vol(fl*(l)) -> J3icf for large- 
order moments (a line in logarithmic scale), where, from a 
numerical fitting, we find that j3i = 0.35 and c\ = 1.9192. 
Therefore, from (0 we have 

E[m k ] x J3i (c\nr) , 

for large k. 

For even-order expected spectral moments (i.e., k - 
s € N), the following holds 



2s for 



E[m 2l .] = -£E[A, 2i ]>-E[A m l x ] 



i=\ 



Define f (n) = n l logn; thus, for any e,5 > (and c\ 
1.9192), we can apply Markov's inequality as follows 

E [^max] 

(c\nr+ erf (n)) 2s 
n E [m 2s ] 



P(X^>(c inr + erf(n)) 2s ) 
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Fig. 3. Comparison between the empirical spectral radius of an RGG 
(circles in the plot) and the values of our analytical upper bound (solid 
line) for n = 1000 and r(n) =d/2n, with expected degrees d= [10:1:100]. 



For large s, one can prove that [8] 



Fornax > cinr + £rf(n)) < «/3iexp sn d log« 

Assuming that s grows as /^n 5 , for j32,<S > 0, we have 



P(Amax >cmr+erf(n)) <«j3iexp 
for all sufficiently large e. Thus, 
. Amax < c\nr 



logn =o(l), 



lim . 



ern 1 5 log« 



L. 



(13) 



In other words, is upper-bounded by cnr + ern i ~ s logn 
with probability 1 for n —> °°. In practice, for a large (but 
finite) «, we can use 1.9192 nr as an upper bound of Amax- 
In Fig. 4, we plot the empirical spectral radius of an RGG 
with n = 1000 and r(n) = d/2n, with expected degrees d = 
[10:1:100] (circles in the figure). We also plot the values of 
our analytical upper bound, 1.9192 nr, in solid line. 

The technique introduced in this subsection is also valid 
for RGG's in higher-dimensions. In general, one can prove 
that for a c/-dimensional RGG that the expected spectral 
moment grows as E {m k \ — > fid (cdn^ . Applying Wigner's 
high-order moment method to this sequence, one can derive 
a probabilistic upper bound similar to ( fT3l . In particular, we 
have that Amax < c^nr' 1 for large n with high probability. In 
the following subsection, we use our results to design the 
connectivity radius of an RGG in order to tame an initial 
viral infection. 

B. Spectral Radius Design 

Once the spectral radius is upper-bounded, our design 
problem becomes trivial. Since (01 represents a sufficient 
condition for local stability around the disease-free state, we 
have the following condition to tame an initial viral infection 
for a rf-dimensional RGG: 

4 ^ 



(c\nr+ erf (n)) 2 



Amax(G(^„;r)) < c d nr" < 



P' 



Fig. 4. Color map representing the evolution of the probabilities of infection 
Pi[n] for = 1,...,1000 in an RGG with n = 1000 nodes, connectivity 
radius r = 0.005, rate of infection ft = 0.020, and recovery rate S = 0.018. 
Each horizontal line represents the value of p, [n] for a particular ;. In this 
color map, blue represents a zero value, green and yellow tones represent 
intermediate values, and red represents values close to one. In this case, we 
observe an epidemic outbreak. 



which implies the following design condition for the con- 
nectivity radius: 

8 y/d 

I , (14) 



J3 c d n 

where c c [ is a positive constant that depends on the dimension 
of T . For example, in the one-dimensional case, we have 
ci = 1.9192; hence, (O becomes r < 5/(1.9192 J3n). We 
now validate this result with several numerical simulations 
of a viral infection in a one-dimensional RGG. 

Consider an RGG with n = 1000 nodes and a connectivity 
radius of r = 0.005 (which implies an average degree of 
10). The resulting spectral radius in this RGG is A ma x = 
17.2629. In our numerical simulations, we choose the initial 
probability of infection to be pi[0] — O.OlUnif [0, 1]; hence, 
approximately 1% of the nodes in the network are initially 
infected. In our first experiment, we choose a rate of infection 
j3 = 0.020, and a recovery rate 8 = 0.018. Since the sufficient 
condition for viral control in (TT~4l > is not satisfied, we cannot 
guarantee an initial infection to be tamed. In Fig. 5 we show 
an image of the evolution of the probability of infection for 
this case. This figure is a color map for the simultaneos 
evolution of pi\n] for / = 1,...,1000. Each horizontal line 
represents the value of p; [n] for a particular ;. In this color 
map, blue represents a zero value, green and yellow tones 
represent intermediate values, and red represents values close 
to one. On the other hand, if we increase the recovery rate to 
S = 0.35 keeping the rest of parameters fixed, we have that 
5//3 = 17.50 > Amax and we satisfy condition OJ. Hence, 
the probability of infection of every node is guaranteed to 
converge towards zero. In Fig. 6, we observe the color map 
for the evolution of the probability of infection in this case, 
where we clearly observe how pi [n] — > for all i. Hence, 
this latter RGG is well-suited to tame initial viral infections. 

V. Conclusions 

In this paper, we have studied the spreading of a viral 
infection in a random geometric graph from a spectral point 



Fig. 5. Color map representing the evolution of the probabilities of infection 
Pi [n] when we increase the recovery rate to S = 0.35 (the rest of parameters 
are the same as we used for Fig. 5). We observe how the probability of 
infection of every node converges towards zero in this case. 



of view. We have focused our attention on studying the 
eigenvalue distribution of the adjacency matrix. We have 
derived, for the first time, explicit expressions for the spectral 
moments of the adjacency matrix as a function of the density 
of nodes and the connectivity radius. We have then applied 
our results to the problem of viral spreading in a network 
with a low-density infection. Using our expressions, we have 
derived upper bounds for the spectral radius of the adjacency 
matrix. Finally, we have applied this upper bound to design 
random geometric graphs that are well-suited to tame an 
initial low-density infection. Our numerical results match our 
predictions with high accuracy. 
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